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ABSTRACT 


Results  from  shallow-water  observational  studies  have  shown  aeoustic  field 
fluetuations  in  excess  of  10  dB  due  primarily  to  non-linear  internal  solitary  waves 
(ISWs).  This  work  concentrates  on  three  limitations  ISWs  pose  to  shallow-water  acoustic 
propagation:  anisotropy,  quasi-deterministic-stochastic  nature,  and  frequency 

dependence.  These  aspects  are  explored  for  low  frequencies  of  75  and  150  Hz  through 
the  development  of  a  full-wave  three-dimensional  parabolic  equation  model  in  which  a 
single  ISW  is  introduced  and  rotated  to  achieve  differing  geometries  between  the  ISW 
and  acoustic  transmission  path.  Two  acoustic  scattering  regimes  were  found  to  exist:  A 
horizontal  refractive  regime  that  generates  very  intense  acoustic  energy  focusing  and 
shadow  zones  near  the  ISW  at  range  when  the  acoustic  transmission  path  is  nearly 
parallel  to  the  ISW  crest,  and  a  vertical  mode  coupling  regime  producing  moderate  to 
strong  acoustic  energy  focusing  and  shadow  zones  for  all  geometries.  These  three- 
dimensional  patterns  are  similar  for  each  frequency,  but  more  intense  for  150  Hz.  The 
results  dramatically  show  ISWs  focus  acoustic  energy,  dependent  upon  frequency  and 
geometry,  which  may  be  exploitable  in  both  a  sonar  performance  and  ambient  noise 
modeling  sense. 
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I.  INTRODUCTION 


Results  from  observational  studies  sueh  as  the  SWARM  (shallow-water  acoustics 
in  a  random  medium)  experiment  in  1995  on  the  continental  shelf  region  of  the  New 
Jersey  coast,  the  New  England  Shelfbreak  PRIMER  from  1995  to  1997,  and  the  Asian 
Seas  International  Acoustics  Experiment  (ASIAEX)  of  2001  in  the  northern  South  China 
Sea  have  illustrated  significant  acoustic  field  fluctuations  in  shallow-water  environments 
with  peak-to-peak  acoustic  energy  fluctuations  of  more  than  10  dB.  An  analysis  of  mode 
coupling  due  to  internal  solitary  waves  (ISWs)  and  their  effects  on  modal  intensity  was 
documented  by  Tielburger,  Einette  and  Wolf  (1997).  Eurther  work  by  Preisig  and  Duda 
(1997)  suggested  ISWs  typically  appear  as  propagating  pulse-like  depressions  of  the 
thermocline  which  cause  time  variations  in  acoustic  intensity  and  phase  attributable  to 
acoustic  normal  mode  coupling.  The  analyses  of  these  observational  studies  demonstrate 
some  of  the  difficulties  in  understanding  shallow-water  acoustic  propagation  and  sonar 
performance  prediction  in  the  presence  of  nonlinear  ISWs. 

ISWs  are  widespread  and  prominent  features  produced  in  shallow-water  regions 
by  tidal  currents  advecting  stratified  water  over  local  bathymetry,  producing  an  internal 
wave  at  tidal  frequency  (an  internal  tide).  The  internal  tide  evolves  non-linearly  into  one 
or  more  oscillations,  appearing  as  a  single  ISW  or  a  packet  of  ISWs.  Many  coastal 
regions  regularly  produce  ISWs  with  amplitudes  of  50  to  100  m  (Jackson,  2004). 

This  thesis  concentrates  on  three  effects  ISWs  impose  on  shallow-water  acoustic 
propagation;  anisotropy,  a  quasi-deterministic-stochastic  nature,  and  frequency 
dependence.  ISWs  impart  anisotropy  to  a  shallow-water  environment  in  that  they  are 
locally  planar  propagating  wave  fronts  containing  strong  horizontal  sound  speed 
gradients  in  the  direction  of  ISW  propagation.  Thus  acoustic  propagation  is  a  function  of 
the  ISW’s  angle  relative  to  the  acoustic  transmission  path,  i.e.,  propagation  paths  along 
verses  across  ISWs  differ  dramatically.  ISWs  are  generally  deterministic,  appearing 
nearly  every  tide  cycle  given  the  same  season,  tidal  phase,  and  bathymetry.  However, 
they  are  stochastic  in  that  their  position,  amplitude,  width,  and  number  vary  based  on 
tidal  phase  and  season.  Thus  the  presence  of  ISW(s)  demonstrates  a  quasi-deterministic- 
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stochastic  shallow-water  acoustic  propagation  environment.  It  is  well  known  that 
shallow-water  aeoustie  propagation  is  strongly  frequency  dependent  owing  to  bottom 
attenuation  and  water  column  absorption.  Due  to  the  small  seale  of  ISWs  and  the 
resonant  interaetion  of  aeoustie  fields  (Young,  2006)  with  ISWs,  there  is  additional 
strong  frequency  dependence  to  the  scattering. 

The  general  influenees  on  shallow-water  aeoustie  propagation  derive  from  the 
seattering  effeets  of  the  surfaee,  bottom,  and  volume.  This  thesis  provides  insight  to 
shallow-water  aeoustie  propagation  concerning  ISW  limitations  by  speoifieally 
addressing  the  dependence  on  volume  scattering  through  development  of  a  full-wave 
three-dimensional  parabolic  equation  (3DPE)  model.  While  volume  scattering  has  been 
modeled  with  techniques  sueh  as  uncoupled  3D  normal  modes,  ray  and  N  x  2D  methods, 
it  has  not  been  treated  with  a  full-wave  3D  modeling  approach. 

The  anisotropie  aspeet  of  ISWs  is  explored  by  observing  the  aeoustie  variability 
as  a  function  of  ISW  orientation.  A  single  ISW  is  introduced  into  the  3DPE  model  which 
is  then  rotated  in  the  xy-plane,  thereby  realizing  differing  angles  between  the  ISW  and 
acoustic  transmission  path.  The  deterministic  character  of  shallow-water  acoustic 
propagation  concerning  ISWs  is  analyzed  by  examining  the  acoustic  energy  field  patterns 
and  energy  flux.  The  stoehastie  character  is  not  explored  in  this  work,  but  can  be  studied 
by  varying  ISW  parameters.  Erequency  dependenee  was  investigated  by  using  75  and  150 
Hz. 

The  results  dramatically  show  that  across-ISW  propagation  is  characterized  by 
moderate  mode  coupling  and  little  out-of-plane  scattering,  or  a  weak  three-dimensional 
effeet,  while  along  ISW  propagation  is  eharaeterized  by  both  strong  mode  eoupling  and 
out-of-plane  scattering,  or  a  strong  three-dimensional  effect  with  clearly  defined 
foeusing/defoeusing  regions.  The  mode  energy  was  found  to  fluctuate  in  excess  of  3  dB 
from  an  environment  free  of  ISWs.  Also,  angular  deflection  of  acoustic  energy  on  the 
order  of  1-5°  was  found  to  exist  near  the  ISW  at  range.  Eurther,  the  angular  defleetion  and 
associated  energy  fluctuations  increase  with  range. 

The  aetual  limitations  caused  by  ISWs  such  as  acoustic  energy  fluctuations  of  10 
dB  and  phase  fluetuations  hampering  signal  proeessing  and  detection  can  be  mitigated. 
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ISWs  shape  the  ambient  noise  field  and  foeus  energy.  Gaining  knowledge  of  these 
patterns  ean  allow  for  deteetion  above  noise.  If  knowledge  is  gained  of  the  phase 
fluetuation  statisties,  better  algorithms  eould  be  developed  to  enhanee  deteetion 
profieieney.  Henee,  this  work  has  implieations  for  both  aetive  and  passive  sonar 
eapabilities  in  addition  to  ambient  noise  modeling. 

Chapter  II  deseribes  the  development  of  the  3DPE  model.  Chapter  III  diseusses 
model  simulations,  data  sets  and  output  display  methods.  Model  analysis  is  provided  for 
aeoustie  transmission  paths  in  strietly  along-  and  aeross-ISW  geometries  in  the  first 
seetion  of  Chapter  IV.  The  second  section  of  Chapter  IV  describes  the  acoustic  energy 
flux  and  resulting  acoustic  modification  regimes.  The  last  section  of  Chapter  IV  provides 
analysis  of  acoustic  transmission  paths  for  intermediate  ISW  geometries.  Conclusions 
and  recommendations  are  provided  in  Chapter  V. 
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II.  DEVELOPMENT  OF  A  THREE-DIMENSIONAL  PARABOLIC 

EQUATION  MODEL 


A,  THREE-DIMENSIONAL  STANDARD  PARABOLIC  EQUATION 

The  standard  parabolic  equation  introduced  by  R.  H.  Hardin  and  F.  D.  Tapper!  in 
the  early  1970s  is  used  to  simulate  acoustic  propagation,  the  derivation  of  which  closely 
follows  that  of  Jensen,  Kuperman,  Porter,  &  Schmidt  (1994).  The  derivation  begins  with 
the  Helmholtz  equation  for  a  homogeneous  density  medium  in  Cartesian  coordinates 

V^p  +  ^p  =  0,  (1) 

c 

where  p  is  the  wave  function  describing  the  acoustic  pressure,  co  the  angular  frequency 
and  c  the  sound  speed.  It  is  assumed  that  the  solution  to  Equation  (1)  takes  the  form 

p{x,y,z)  =  i//{x,y,z)e‘'‘'‘\  (2) 

where  ip  is  the  complex  wave  amplitude  and  =  (u/cg  a  reference  wavenumber. 

Taking  the  first  and  second  derivatives  of  the  solution  yields 

p'  =  (3) 

and 

p"  =  .  (4) 

Substituting  Equations  (3)  and  (4)  into  the  Helmholtz  equation  yields 


d^y/  d^y/  d^y/ 
dx^  dy^  dz^ 


+  lik^y/'  -kly/  + 


y/  =  0 . 


(5) 


The  wave  is  assumed  to  be  slowly  changing  in  x.  Hence,  the  paraxial  approximation 
expressed  by 


d^yr 


« lik^y/ 


(6) 


is  then  introduced  to  obtain  the  form 
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(V) 


likr. 


dx  dy  dz 


The  term  -  k^  can  be  manipulated  by  using  the  index  of  refraction 


n  = 


(8) 


where 


l-% 

V  c  y 


[kl-k^)  =  kl  1 — ^  =-kl[n-\) 


(9) 


to  yield  the  final  form  of  the  three-dimensional  standard  parabolic  equation 


dy/  d\  0V  ,  ,2/  2 


U  ^  ^2  ^2 

ox  qy  oz 


-  + 


k\  {ji^  -  \^y/  =  0 . 


(10) 


We  reiterate  here  that  the  numerical  modeling  of  this  shallow-water  problem  treats  the 
seabed  as  a  fluid  medium,  thus  neglecting  density  changes  in  depth. 


B,  MODELING  METHOD 

Range  independent  models  treat  the  environment  in  a  homogeneous  sense,  where 
range  dependent  models  are  capable  of  handling  environmental  variations  in  range.  The 
technique  used  in  this  work’s  model  involves  using  a  three-dimensional  parabolic 
equation  (3DPE)  to  simulate  the  pressure  field  through  a  range  dependent  water  column. 
The  3DPE  model  as  developed  does  not  account  for  spatial  bottom  variability,  temporal 
and  spatial  variability  of  the  surface,  or  density  variability. 

Described  by  Knepp  (1983),  a  multiple  phase  screen  propagation  technique  is 
employed,  whereby  a  given  volume  is  divided  into  a  number  of  thin  layers,  each  being 
perpendicular  to  the  direction  of  propagation.  The  irregularities  in  each  layer  are  replaced 
by  the  geometric  optics  phase  change  imparted  after  propagation  through  the  layer,  or  the 
phase  screen. 

The  split-step  Eourier  algorithm  technique  has  been  utilized  extensively  to  solve 

the  PE  and  is  likewise  used  to  solve  the  3DPE  in  this  work.  Jensen  et  al.  (1994)  derives 
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the  split-step  marching  solution  in  cylindrical  coordinates.  This  work  uses  a  Cartesian 
coordinate  form  such  that 


y/{x,y,z)  =  F^-\ 


iAx  ^ 


F 


w(= 


0’ 


3^0  5  ^0 


(11) 


where  a  Fourier  transform  in  y  and  z  is  applied  to  the  original  wave  function  and  phase 
screen  functions,  y/[xQ,y^,ZQ^  and  expj— (xo,yQ,Zo)-lJ Axj  respectively.  The 


final  result  is  obtained  by  applying  the  inverse  Fourier  transform  to  the  aforementioned 

calculation  and  free  propagator  function  exp  \  U  . 

Ikr, 


The  stability  of  this  method  concerning  a  2DPE  model  was  confirmed  for  75,  150 
and  300  Hertz  by  Young  (2006).  Results  from  the  3DPE  model  were  then  compared  to 
those  of  the  2D  case  to  ensure  accuracy.  Due  to  computational  restrictions,  this  work 
provides  an  analysis  of  75  and  150  Hz  only. 

Proper  computational  parameters  for  the  3DPE  model  were  then  established  as 
seen  in  Table  1.  Of  importance  were  the  lateral  and  depth  grid  spacings,  (r/y)  and  (Jz) 

respectively,  and  two  march-step  intervals  (r&)  and  (dhx) .  If  parameter  values  were  too 

large,  the  simulated  acoustic  pressure  field  lacked  accuracy.  If  too  small,  computing  time 
and  memory  constraints  became  an  issue. 

Two  march-step  intervals  were  used  in  the  3DPE  model  to  optimize  computing 
time.  The  first,  (r/x)  ,  was  chosen  as  the  base  march-step  interval  and  was  found  to 

provide  accurate  results  when  compared  to  a  similar  2DPE  model.  The  second,  (r&x) , 

smaller  march-step  was  established  primarily  to  provide  a  smaller  interval  over  which 
energy  flux  intensities  could  be  evaluated.  The  smaller  march-step  was  applied  once  for 
every  base  march-step.  This  considerably  decreased  computing  time  as  compared  to 
using  the  smaller  march-step  exclusively.  The  computational  parameters  chosen  for  the 
3DPE  model  in  this  work  are  provided  in  Table  1. 
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Base  march-step  interval 

dx 

10  m 

Small  march-step  interval 

dxx 

1  m 

Lateral  grid  spacing 

dy 

0.73  m 

Depth  grid  spacing 

dz 

0.73  m 

Range 

X 

5000  m 

Lateral  extent 

ymax 

3000  m 

Maximum  computational  depth 

D 

375  m 

Table  1 .  Grid  spacing  and  march-step  interval  parameters. 


C.  BOUNDARY  CONDITIONS 

The  ocean  surface  was  modeled  as  a  pressure  release  boundary  using  an  image 
ocean.  This  imposed  the  constraint  of  zero  pressure  at  the  surface  y/  [x,  y,  O)  =  0 .  The 

bottom  and  lateral  boundary  conditions  were  absorbing  in  nature.  Boundary  attenuation 
was  achieved  by  introducing  artificial  sponge  layers  in  such  a  way  that  all  waves  were 
totally  absorbed  and  no  reflection  occurred.  This  was  accomplished  by  imposing  a 
gradual  loss  of  amplitude  as  y  approached  the  lateral  limits  and  z  approached  the  bottom. 
The  form  of  the  loss  at  each  step  is 


L(y,z)  =  exp< 

-(ddx 

^  V  '^.kmax  J  ^  ^  ' 

J, 

where  a  =  0.08,  (5  =  0.04  and  D  is  the  maximum  computational  depth.  This  form 
effectively  prevents  any  acoustic  energy  from  penetrating  ay.D  -  above  the  bottom 
and  likewise,  a  x  y^^^  -12  m  from  the  lateral  boundaries. 


D,  SOURCE/STARTING  FIELDS 

The  3DPE  model  uses  an  analytical  source  function  designed  to  match  the  farfield 
result  for  a  point  source  solution  of  the  Helmholtz  equation  in  a  homogeneous  medium. 
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The  result  of  the  derivation  of  a  Gaussian  point  souree  found  in  Jensen  et  al.  (1994)  was 
used  in  this  work: 


y/[x  =  Q,y  =  Q,z)  =  k^e  ^  ,  (13) 

where  is  the  effeetive  souree  level  and  is  the  souree  depth.  The  value  of  k^  is 
obtained  by 


(14) 


where  /g  is  the  frequeney  and  Cg  is  the  referenee  sound  speed.  In  this  work,  the  referenee 
sound  speed  was  1500  m/s  and  the  souree  depth  50  m. 


E.  EXTRACTION  OF  MODE  AMPLITUDES 

Normal  mode  amplitudes  were  ealculated  for  modes  one  through  twenty  in  this 
work.  They  were  extraeted  by  expanding  the  funetion  ^(x,y,z)  in  terms  of  the 

unperturbed  modes,  ,  yielding 


y/{x,y,z)  =  ^A,Xx,y)(l)„{z) 


(15) 


n=\ 


and 


=0, 


(16) 


where  {x,y^  are  the  mode  amplitudes  and  k  =  cojc  (z)  (Colosi,  2006). 


The  mode  amplitude  equation  is  obtained  by  substituting  equation  (15)  into  the 
three-dimensional  parabolic  equation  (10),  multiplying  by  the  mode  function  (p^  and 
integrating  over  depth  to  yield 


^(Pn  I  i 
dx  2k^  dy^ 


N 


m=l 


(17) 


where 
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(18) 


Pmn  y)=K  \  , 

-D  ^0 

ft  (•■'•->’)  =  4  (■■'.J')'?""''. 

K=KI^K. 


mn  m  n 


(19) 

(20) 

(21) 


and  D  is  the  computational  depth. 

Rather  than  attempt  the  solution  for  Equation  (17),  the  mode  amplitudes  were 
obtained  by  solving  Equation  (10)  for  the  pressure  field  ^(x,y,z) ,  and  multiplying  by 

the  unperturbed  modes,  (j)^ ,  then  depth  integrated  to  yield  (Colosi,  1996) 


D 

A{x,y)  =  \v^{x,y,z)^^{z)dz . 
0 


(22) 


F.  ENERGY  FLUX  INTENSITY 

The  energy  flux  intensity  was  calculated  for  modes  one  through  four  in  this  work. 
Frisk  (1994)  defines  the  energy  flux  intensity  as 

(23) 


where  is  a  reference  density.  For  normal  modes  the  vertical  flux  is  effectively  zero,  or 
1^=0.  Thus  the  only  concern  is  with  components  in  the x  andy  directions.  For  notational 
efficiency  it  is  useful  to  define 


(24) 


Expanding  equation  (23)  in  terms  of  the  unperturbed  modes  yields 


/  = 


f  N  \  f  N  \  f  N  \  f  N  \ 

yas  V  -  y^v  V  yas 

nr  n  mr  m  mrm  nr  n 

\  n=\  )  \m=\  J  Vw=l  J  V  w=l  J 


(25) 
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Equation  (25)  now  yields 


I  =■ 


(  N  ^  f  ^  ^  ^ 

-  E4>.  Ev.4A 

V  «=1  J\m=\  J  Vw=l  J\n=\ 


which  when  depth  integrated  becomes 


I  = 


Acop^ 


EE  4  j  -f  E  E  4: 

«=1  m=\  0  J  y  m=l  n=\  0  J 


Because  of  mode  orthogonality,  this  results  in 

^fN  \  f  N 


I  = 


4®A 


E4V,4‘  -  E4‘v,4 

V  «  J  \  n 


(26) 


(27) 


(28) 


leading  to  the  modal  energy  flux  intensity  equation  used  in  the  3DPE  model: 


in{^^y)  = 


A(op^ 


Ay, A -Ay, A 


(29) 


Note  that  the  small  range  step  dxx  discussed  in  section  B  is  used  to  numerically  compute 
the  range  derivatives  d/dx  in  Equation  (29). 


G.  SOUND  SPEED  FIELD 

The  environment  used  in  this  work  was  a  range  independent  background 
waveguide  consisting  of  four  layers:  Surface  layer,  thermocline,  below  thermocline 
isovelocity  layer  and  an  isovelocity  bottom  layer.  Transitions  between  these  layers  were 
treated  with  a  hyperbolic  tangent.  The  unperturbed  sound  speed  field  was  obtained  using 
the  following  equation: 


r 

c,  - 

C  —  C; 

5  / 

1  +  tanh 

^-^th 

\ 

+ 

Ci-Ck 

( 

1  +  tanh 

z-(z,+4E)T 

s 

V 

2  . 

-  A,  _ 

) 

2 

V 

L  ^  JJ 

where  is  the  surface  sound  speed,  c,  the  below  thermocline  sound  speed,  c* 
the  bottom  sound  speed,  the  bottom  depth,  the  central  depth  of  the  thermocline, 
the  thermocline  width,  and  L  the  transition  length  scale.  The  transition  length  scale 
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is  needed  beeause  of  the  Fourier  split-step  numerieal  method.  Figure  1  depiets  the 
unperturbed  sound  speed  profde.  Table  2  displays  the  values  used  in  the  3DPE  model. 


Surfaee  sound  speed 

G 

1520  m/s 

Below  thermoeline  sound  speed 

G 

1480  m/s 

Bottom  sound  speed 

G 

1980  m/s 

Bottom  depth 

75  m 

Central  thermoeline  depth 

G/, 

30  m 

Thermoeline  width 

A,, 

15  m 

Transition  length  seale  75  Hz 

^15Hz 

20  m 

Transition  length  seale  150  Hz 

10m 

Table  2.  Sound  speed  field  parameters. 


H.  INTERNAL  SOLITARY  WAVE 

Work  by  Young  (2006)  has  held  that  using  a  Gaussian  form  of  a  soliton  ean 
aeeurately  approximate  the  hyperbolie-seeant-squared  form  of  a  soliton  whieh  is  an  exaet 
solution  to  the  Kortweg  de  Vries  wave  equation.  Following  Young  (2006)  we  introduee  a 
single  soliton,  or  ISW  to  the  sound  speed  field  using  the  Gaussian  form 


r 

4'(x,y,z;Xo,yo,6»)  =  ^sin 


j 


1  -  tanh 


L 


■exp 


[(x-Xo)sin6'  +  (y-yo)eos6»]^ 

a! 


(31) 

where  is  the  maximum  soliton  amplitude,  z  is  depth,  z^  the  bottom  depth,  L  the 
transition  length  seale,  the  soliton  width,  6*  the  angle  between  the  x-axis  and  ISW, 
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and  (xq,  Jq)  the  hinge  point  about  whieh  the  ISW  rotated.  While  the  3DPE  method  eould 

easily  aeeommodate  the  proper  hyperbolie-seeant-squared  soliton  shape,  we  used  the 
approximate  Gaussian  shape  so  that  eomparisons  eould  be  made  to  the  mode  eoupling 
theory  in  Young  (2006).  The  perturbed  sound  speed  field  eontaining  the  ISW  is  then 
obtained  by 

c(x,y,z)  =  c(z)  +  —  ^(x,y,z;xo,yo,0) .  (32) 

oz 

Figure  2  displays  the  sound  speed  field  in  the  xz-plane  containing  a  single  ISW  oriented 
parallel  to  they-axis  at  x  =  1000  meters. 


I.  MODEL  VALIDATION 

The  3DPE  model  was  validated  in  three  phases.  The  first  test  consisted  of  using  a 
point  source  to  evaluate  the  spherical  spreading  along  the  x-axis.  The  spherical  spreading 
was  found  to  be  satisfactory.  The  second  test  consisted  of  using  a  point-mode  source  to 
evaluate  cylindrical  spreading  and  the  conservation  of  mode  energy.  The  cylindrical 
spreading  was  found  to  be  adequate  and  numerical  coupling  was  nonexistent  for  the 
homogeneous  case.  In  the  third  test,  the  results  for  an  ISW  angle  of  90°  (directly  across 
ISW  transmission  path)  were  compared  to  that  of  Young’s  2D  calculations  using  similar 
soliton  parameters  and  good  agreement  was  found. 
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III.  MODEL  SIMULATIONS  AND  DISPLAYS 


A.  MODEL  SIMULATIONS 

Twenty  separate  3DPE  model  simulations  were  completed  for  each  frequency,  75 
and  150  Hz.  The  first  3DPE  simulation  was  completed  for  an  unperturbed  environment 
with  no  ISW  referred  to  as  the  homogeneous  case.  Rotational  limiting  cases  were  then 
defined  and  simulated.  The  lower  limiting  case  was  defined  as  the  0“  ISW  angle  at  which 
the  ISW  was  oriented  along  the  x-axis.  The  upper  limiting  case  was  defined  as  the  90° 
ISW  angle  at  which  the  ISW  was  oriented  perpendicular  to  the  x-axis  and  centered  at 
X  =  1000  m  .  Easily,  3DPE  model  simulations  were  completed  every  5°  between  the  two 
limiting  cases.  Eigure  3  displays  the  ISW  locations  used  in  the  model  simulations. 

Eor  each  simulation  the  mode  amplitudes  were  computed  as  well  as  the 

energy  flux  vectors  •  For  the  purposes  of  displaying  the  change  in  mode  energy 

for  the  various  cases  we  formed  the  ratio  (x,y,6’)p/|H„Q  (x,y)|^ ,  where  |H„o(x,y)|^ 

was  the  mode  energy  in  the  homogeneous  case.  In  this  analysis  we  focus  on  modes  1-4, 
although  mode  4  of  the  150  Hz  data  sets  were  found  to  contain  very  little  energy  and 
became  unreliable  when  normalized  by  the  homogeneous  energy. 


B,  MODEL  DISPLAYS 

The  model  results  in  this  work  are  presented  by  two  methods;  5250  m  x  3000  m 
plan  view  plots,  as  in  Eigure  3,  displaying  the  azimuthal  dependence  of  normal  mode 
energy  fluctuation  patterns  (focusing/defocusing)  and  refractive  effects,  and  waterfall 
plots  that  demonstrate  the  same  three-dimensional  effect  as  a  function  of  ISW  angle  for  a 

given  range.  In  the  second  method,  displays  were  generated  for  (x, 7,6*0)  for 
X  =  2,  2.5,  3,  3.5,  4,  4.5  and  5  km.  All  plan  view  plots  denote  the  position  of  the  ISW 
with  a  black  line.  All  plots  were  standardized  to  display  the  same  intensity  range  for 
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IV.  MODEL  ANALYSIS 


A.  ACOUSTIC  SCATTERING  PHYSICS  FOR  LIMITING  CASE 

SIMULATIONS 

1,  75  Hz  Normal  Mode  Energy  Interactions  for  ISW  Limiting  Cases 

Figure  4  depicts  75  Hz  modal  energy  patterns  for  the  case  in  which  the  ISW  is 
propagating  in  the  y-direction  so  that  the  crest  of  the  ISW  is  entirely  along  the  x-axis 
defining  the  0°  ISW  angle  limiting  case.  In  this  orientation  there  is  no  gradient  of  sound 
speed  in  the  x-direction  and  the  gradient  is  strongest  in  the  y-direction;  clearly,  this  case 
only  allows  refraction  and  mode  coupling  in  the  y-direction.  In  these  figures  a  clear 
negative  lensing  effect  is  seen  such  that  along  the  ISW  axis  (x-axis)  there  is  a  significant 
and  widening  shadow  zone  of  low  intensity.  Also,  there  is  a  strong  dependence  on  mode 
number  for  this  diverging  lens  pattern  with  the  largest  change  in  modal  energy  seen  for 
mode  2.  Figure  5  shows  the  fractional  change  in  mode  energy  along  the  shadow  zone  axis 
(x-axis)  and  it  is  clear  that  significant  mode  coupling  is  occurring,  especially  between 
modes  1  and  2  where  the  mode  energy  appears  to  oscillate  back  and  forth  between  the 
modes.  Also  note  that  the  spherical  spreading  of  the  mode  energies  has  been  taken  out  of 

this  display  because  we  normalize  by  •  Thus  the  decrease  in  energy  with  range 

seen  in  Figure  5  is  due  to  the  negative  gradients  in  sound  speed  experienced 
perpendicular  to,  and  on  either  side  of  the  ISW.  The  acoustic  energy  flux,  covered  in  part 
B  of  this  analysis,  indicates  refraction  angles  up  to  5°  near  the  ISW.  Thus,  in  this 
geometry  there  is  a  clear  combination  of  scattering  mechanisms  involving  both  refraction 
and  mode  coupling. 

Figure  6  displays  the  modal  energy  patterns  for  the  case  in  which  the  ISW  is 
propagating  in  the  x-direction,  so  that  the  crest  is  entirely  parallel  to  the  y-axis  defining 
the  90°  ISW  angle  limiting  case.  In  this  orientation  there  is  no  gradient  of  sound  speed  in 
the  y-direction  and  the  gradient  is  strongest  in  the  x-direction;  clearly  this  case  only 
allows  mode  coupling  in  the  x-direction.  The  acoustic  energy  flux  shows  a  complete 
absence  of  refraction  in  this  geometry.  Thus,  there  is  no  3D  effect  exhibited  in  these 

figures  since  the  mode  energy  simply  shifts  as  the  acoustic  field  passes  the  1000m  range. 
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Figure  7  shows  the  fraetional  ehange  in  mode  energy  along  the  x-axis  demonstrating 
moderate  mode  eoupling,  partieularly  between  modes  2  and  3. 

2,  150  Hertz  Normal  Mode  Energy  Interactions  with  ISW  for  Limiting 

Cases 

For  the  150  Hz  simulations,  Figures  8  and  9  show  signifieant  initial  modal  energy 
in  the  first  three  modes,  while  mode  4  is  very  weak  (see  espeeially  Figure  9).  In  eontrast, 
the  first  four  modes  contained  significant  mode  energy  in  the  75  Hz  simulations  as  seen 
in  Figures  5  and  7.  As  in  the  75  Hz  simulations,  strong  mode  coupling  exists  between  the 
first  three  modes  with  gradual  energy  depletion  in  all  modes  along  the  axis  of  the  ISW  in 
the  along  ISW  transmission  path  with  moderate  mode  coupling  evident  between  modes 
two  and  three  after  interaction  with  the  ISW  in  the  strictly  across-ISW  transmission  path. 
Figure  10  shows  that  the  symmetric  3D  effect  for  the  150  Hz  case  in  the  along  ISW 
acoustic  transmission  path  is  similar  to  that  of  the  75  Hz  case. 


B,  ACOUSTIC  ENERGY  FLUX 

Although  three-dimensional  effects  are  observed  concerning  all  modes,  they  are 
most  dramatic  and  best  expressed  by  mode  2.  Hence,  figures  representing  characteristic 
behavior  are  provided  for  mode  2  exclusively  in  the  remainder  of  this  work. 

As  the  ISW  is  rotated,  two  distinct  regimes  become  apparent  based  on  an 
evaluation  of  the  acoustic  energy  flux.  Figure  11  displays  the  75  Hz  mode  2  energy 
fluctuation  for  an  ISW  angle  of  5°.  The  boxed  regions  refer  to  Figures  12  and  13  which 
illustrate  the  acoustic  energy  flux  vectors  for  both  the  unperturbed  (no  ISW)  simulation 
and  the  perturbed  (ISW)  case.  Figure  12  shows  only  a  change  in  magnitude  between  the 
unperturbed  and  perturbed  acoustic  fluxes  indicating  mode  coupling  is  the  only 
mechanism  responsible  for  focusing  modal  energy  through  the  ISW.  In  contrast.  Figure 
1 3  displays  both  magnitude  and  directional  changes  in  the  acoustic  fluxes  indicating  that 
refraction  and  perhaps  mode  coupling  are  leading  to  the  intense  focusing  and  defocusing 
along  the  ISW. 
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Figure  14  displays  the  energy  fluetuation  and  refraetion  angles  ealeulated  from 
the  differenees  in  aeoustie  flux  between  the  unperturbed  and  perturbed  simulations  for  75 
and  150  Hz  simulations  with  a  5°  ISW  angle.  The  strongest  gradient  in  sound  speed  lies 
near,  and  is  perpendieular  with,  the  ISW.  Thus,  the  extent  of  the  refraetive  effeet  is 
limited  to  regions  near  the  ISW  at  range.  Further,  the  highest  sound  speed  is  along  the 
erest  of  the  ISW  eausing  negative  sound  speed  gradients  on  either  side  of  the  erest.  Thus, 
the  resultant  refraetion  angles  are  negative  in  the  -y  direetion  from  the  ISW  and  positive 
in  the  +y  direetion.  Further,  refraetion  of  modal  energy  oeeurs  only  for  ISW  angles  of  less 
than  50°.  Figure  15  illustrates  the  last  vestige  of  the  refraetive  effeet  at  an  ISW  angle  of 
45°  in  the  75  Hz  simulation,  while  the  refraetive  effeet  is  negligible  eoneerning  the  150 
Hz  ease.  This  diseovery  that  refraetion  is  important  up  to  a  45°  ISW  angle  is  new  and 
mueh  larger  than  previously  predieted  (Lyneh,  Colosi,  Gawarkiewiez,  Duda,  Pieree, 
Badiey,  Katsnelson,  Miller,  Siegmann,  Chiu  &  Newhall,  2006). 

Refraetion  of  mode  energy  oeeurs  only  near  the  ISW  and  for  lower  ISW  angles 
whieh  eorresponds  to  along  ISW  aeoustie  transmission  paths.  The  refraetion  does  not 
appear  to  be  frequeney  dependent.  However,  given  that  only  two  low  frequeney 
simulations  were  eompleted  a  more  eomplete  analysis  is  warranted  to  justify  this 
statement.  Mode  eoupling  appears  to  oeeur  at  all  ISW  angles  eorresponding  to  all 
aeoustie  transmission  paths,  although  the  eoupling  is  more  intense  for  those  paths  along 
the  ISW.  Further,  the  mode  eoupling  behavior  is  frequeney  dependent. 


C.  ACOUSTIC  SCATTERING  PHYSICS  FOR  ROTATIONAL 

SIMULATIONS 

1.  75  Hz  Normal  Mode  Energy  Interactions  with  ISW  for  Rotational 

Cases 

Figure  4  depiets  strong,  symmetrie  mode  energy  foeusing  behavior  due  to 
refraetion  and  mode  eoupling  for  aeoustie  transmission  paths  along  the  ISW.  As  the  ISW 
angle  inereases,  the  foeusing  in  the  refraetive  regime  beeomes  more  intense  as  seen  in 
Figure  16.  The  very  intense  foeusing  along  the  souree  side  of  the  ISW  is  aeeompanied  by 
a  strong  shadow  region  on  the  non-souree  side  of  the  ISW.  The  upper  right  panel  of 
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Figure  16  depicts  the  mode  energy  fluctuation  with  a  15“  ISW  angle  and  is  representative 
of  the  gradual  constriction  of  the  source  side  focusing  region.  The  lower  left  panel  of 
Figure  16  displaying  the  mode  energy  fluctuation  at  an  ISW  angle  of  35“  illustrates  the 
continuing  constriction  of  this  focusing  region  with  increasing  ISW  angle.  The  along- 
ISW  focusing  region  becomes  virtually  nonexistent  in  response  to  the  decrease  in  the 
refractive  effect  as  seen  in  the  lower  right  panel  of  Figure  16  displaying  the  mode  energy 
fluctuation  at  an  ISW  angle  of  65“. 

Figure  16  also  displays  regions  of  alternating  focusing  and  defocusing  on  the  non¬ 
source  side  of  the  ISW.  This  effect  is  caused  by  mode  coupling  dependence  on  the 
relative  phase  of  modes  when  they  encounter  the  soliton  (Young,  2006). 

2,  150  Hertz  Normal  Mode  Energy  Interactions  with  ISW  for  Rotational 

Cases 

Figures  16  and  17  illustrate  the  progression  of  the  3D  acoustic  energy  pattern  as 
the  ISW  angle  increases  from  5“  through  65“  for  75  and  150  Hz,  respectively.  Clearly,  the 
energy  fluctuation  patterns  of  the  75  and  150  Hz  cases  are  similar;  however  the  effect  of 
the  ISW  on  acoustic  energy  focusing  is  stronger  at  150  Hz.  The  magnitude  of  the  focused 
150  Hz  acoustic  energy  in  the  refractive  regime  along  the  ISW  is  25%  higher  than  that  of 
the  75  Hz  case.  Further,  this  frequency  dependence  is  not  confined  to  the  refractive 
regime,  but  is  manifest  throughout  the  simulated  waveguide,  with  the  same  25%  increase 
in  focused  acoustic  energy  found  in  the  mode  coupling  regime  on  the  non-source  side  of 
the  ISW. 


3,  Energy  Focusing/Defocusing  as  a  Function  of  ISW  Angle 

To  gain  insight  as  to  focusing/defocusing  of  normal  mode  energy  given  an 
acoustic  source  and  ISW  angle,  it  is  useful  to  depict  the  normal  mode  energy  fluctuation 
as  a  function  of  the  ISW  angle  as  explained  in  Chapter  III.  Figure  18  illustrates  the  75  Hz 
mode  2  energy  fluctuations  at  ranges  of  2.5,  3.5  and  4.5  km. 

The  greatest  effect  occurs  for  ISW  angles  of  less  than  20“.  However,  there  is  a 
marked  effect  near  the  ISW  for  ISW  angles  up  to  45“  and  to  a  lesser  extent  for  the 
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remaining  ISW  angles  up  to  the  90°  limiting  ease.  Clearly  illustrated  is  the  appreeiable 
effeet  of  the  refraetive  regime  for  ISW  angles  of  less  than  20°  noted  by  the  high 
variability  of  intensity  in  the  upper  left  eorner  of  eaeh  figure. 

Figure  19  displays  the  150  Hz  mode  2  energy  fluetuations  at  similar  range  for 
eomparison.  The  eongrueney  of  the  pattern  and  the  ISW’s  effeet  in  the  refraetive  regime 
are  very  similar.  However,  in  stark  eontrast  to  the  75  Hz  ease  is  the  inerease  in  intensity 
within  the  mode  eoupling  regime  seen  in  the  upper  right  of  eaeh  figure.  While  the  mode 
eoupling  effeet  was  appreeiable  at  75  Hz  at  mid-ISW  angles  (20°  to  60°),  the  effeet  is 
mueh  more  dramatie  at  150  Hz. 
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V.  CONCLUSIONS 


Although  this  work  considered  only  a  single  ISW  with  varying  geometries  and 
two  diserete  frequeneies  in  one  shallow-water  aeoustie  waveguide,  important  links 
between  the  anisotropy  imparted  by  ISWs  and  the  resulting  frequeney  dependent  3D 
effeets  have  been  made.  This  work  has  established  that  two  aeoustie  seattering  regimes 
exist.  The  refraetive  regime  generates  very  intense  aeoustie  energy  foeusing  and  shadow 
zones  near  the  ISW  but  only  for  geometries  where  the  aeoustie  transmission  path  is 
nearly  along  the  ISW  erest  (0°  to  20°).  In  eontrast,  the  mode  eoupling  regime  produees 
strong  aeoustie  energy  foeusing  and  shadow  zones  for  nearly  all  geometries,  but 
partieularly  for  along  ISW  aeoustie  transmission  paths. 

Modeling  teehniques  earried  out  only  in  separate  vertieal  planes  where  3D  and 
mode  eoupling  effeets  are  negleeted  may  be  aeeurate  for  strietly  aeross  ISW  propagation 
paths.  However,  this  work  has  shown  that  even  a  small  departure  from  this  geometry 
results  in  inereased  mode  eoupling  and  3D  effeets.  Deterministieally,  the  3DPE  model 
reveals  intensifieation  patterns  that  are  only  implied  by  2D  or  A  x  2D  methods. 

The  ISW’s  effeet  on  the  overall  3D  aeoustie  energy  pattern  was  not  found  to  be 
frequeney  dependent.  However,  the  ISW’s  effeet  on  the  magnitude  of  the  foeusing  is 
frequeney  dependent.  To  wit,  the  150  Hz  simulations  display  a  25%  inerease  in  the 
magnitude  of  the  foeused  energy  over  simulations  of  the  lower  75  Hz  frequeney.  In  a 
relative  sense,  the  effeet  of  this  behavior  is  realized  at  smaller  angular  departures  from 
the  aeross-ISW  aeoustie  transmission  path  for  150  Hz. 

Normal  mode  eoupling  was  shown  to  go  from  lower-order  modes  to  higher-order 
modes  and  viee-versa.  The  analysis  demonstrates  eoupling  behavior  and  intensifieation 
patterns  are  governed  by  the  geometry  of  the  aeoustie  transmission  path  and  ISW.  In  a 
praetieal  sense,  ISWs  foeus  aeoustie  energy  dependent  upon  frequeney  and  geometry, 
whieh  may  be  exploitable  in  both  a  sonar  performanee  and  ambient  noise  modeling  sense. 

Realistieally,  3D  intensifieation  and  mode  eoupling  effeets  oeeur  eontinually  as 
an  aeoustie  field  interaets  with  environmental  inhomogeneities.  This  work  has  shown  the 
eomplexhy  of  shallow-water  aeoustie  propagation  in  the  presenee  of  a  single  ISW. 
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Future  work  should  explore  the  stochastie  nature  of  anisotropy  within  the  acoustic 
waveguide  by  varying  the  number,  shape  and  amplitude  of  ISWs.  Further,  establishing  a 
realistic  acoustic  waveguide  by  including  bathymetry  and  resolving  proper  surface  and 
bottom  attenuation  parameters  could  provide  practical  results  for  a  comprehensive 
comparison  to  observational  studies  such  as  the  shallow-water  acoustic  experiment 
(SW06)  off  the  New  Jersey  coast  in  late  summer  2006. 
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FIGURES 


Unperturbed  Sound  Speed  Profile 


Figure  1 .  Background  sound  speed  profde,  c  (z) ,  used  in  the  3DPE  model. 
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Sound  Speed  (m/s)  vs  Range  and  Depth 


Range  (m) 

Figure  2.  Perturbed  sound  speed  field  in  the  xz-plane  centered  at  j  =  0.  A  single 

ISW  is  located  at  x  =  1000  m  with  an  amplitude  of  10  m,  width  of  100  m,  and  an  angle  of 
90°. 
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ISW  Angles  Computed 


Range  (m) 


Figure  3.  ISW  orientations  used  in  the  3DPE  model.  Model  simulations  were 

eompleted  every  five  degrees  for  ISW  angles  of  zero  to  ninety  degrees. 
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Figure  4.  Normal  mode  energy  fluetuation,  |4,|  / |4io|  >  FIz  for  ISW  angle  of 

0°  and  modes  1-4.  Mode  1  is  shown  in  the  upper  left  panel,  mode  2  upper  right,  mode  3 
lower  left,  and  mode  4  lower  right. 
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Normal  Mode  Energy  (|a^P) 


75  Hz  Normal  Mode  Energy  0°  ISW  Angle 


Figure  5.  Normal  mode  energy,  (x,^  =  O)  ,  at  75  Hz  for  0°  ISW  angle  and 

modes  1-4. 
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Figure  6.  Normal  mode  energy  fluetuation,  |4,|  / |4,o|  ’  75  Hz  for  ISW  angle  of 

90°  and  modes  1-4.  Mode  1  is  shown  in  the  upper  left  panel,  mode  2  upper  right,  mode  3 
lower  left,  and  mode  4  lower  right. 
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Normal  Mode  Energy  (|a^P) 


75  Hz  Normal  Mode  Energy  90°  Grazing  Angle 


Figure  7.  Normal  mode  energy,  (x,^  =  O)  ,  at  75  Hz  for  90°  ISW  angle  and 

modes  1-4. 
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Figure  8. 
modes  1- 


150  Hz  Normal  Mode  Energy  0°  ISW  Angle 


Normal  mode  energy,  (x,^  =  O)  ,  at  150  Flz  for  0°  ISW  angle  and 

4. 
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Normal  Mode  Energy  (|a^P) 


150  Hz  Normal  Mode  Energy  90°  ISW  Angle 


Figure  9.  Normal  mode  energy,  (x,y  =  O) ,  at  75  Hz  for  90°  ISW  angle  and 

modes  1-4. 
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Cross  Rangs  (m)  Cross  Range  (m) 


19D  Hz  Model  Energy  Fluctuation  0°ISW  Angle 


Range  (m) 

193  Hz  Mode  3  Energy  Fluctuation  0°ISW  Angle  (13^1^/13^1^) 


Range  (m) 


150  Hz  Mode  2  Energy  Fluctuation  0°  ISW  Angle  (la^^p/ia^p) 


Range  (m) 
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Figure  10.  Normal  mode  energy  fluctuation,  |4,|  / |4,o|  >  at  150  Hz  for  ISW  angle  of 

0°  and  modes  1-4.  Mode  1  is  shown  in  the  upper  left  panel,  mode  2  upper  right,  and  mode 
3  lower  left. 
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75  Hz  Mode  2  Energy  Fluctuation  5°  ISW  Angle  (|a^p/|a^p) 


Range  (m) 

Figure  11.  Normal  mode  energy  fluctuation,  |4,|  / |4io|  >  Flz  for  ISW  angle  of 

5°  and  mode  2.  The  boxed  regions  correspond  to  Figures  12  and  13. 
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75  Hz  Mode  2  Energy  Fluctuation  5°ISWAngle  (|a^p/|a^p) 
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Figure  12.  Normal  mode  energy  fluctuation,  |4,|  / |4io|  >  Flz  for  ISW  angle  of 

5°  and  mode  2,  overlaid  with  the  perturbed  and  unperturbed  acoustic  energy  flux  vectors 
illustrating  a  mode  coupling  regime. 
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75  Hz  Mode  2  Energy  Fluctuation  5°  ISW  Angle  (|a^p/|a^p) 
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Figure  13.  Normal  mode  energy  fluctuation  ,  at  75  Hz  for  ISW  angle  of 

5°  and  mode  2,  overlaid  with  the  perturbed  and  unperturbed  acoustic  energy  flux  vectors 
illustrating  a  refractive  regime. 


150  Hz  5°  ISW  Angle  Mode  2  Energy  Fluctuation  (|ap,r/|aj.|r)  Refraction  Angle  (deg) 
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Figure  14.  Normal  mode  energy  fluctuation,  |4,|  /lAol  >  for  ISW  angle  of  5°  and 

mode  2,  overlaid  with  the  refraction  angle  in  degrees.  The  75  Hz  case  is  shown  in  the  left 
panel  and  150  Hz  on  the  right. 
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150  Hz  45°  ISW  Angle  Mode  2  Energy  Fluctuation  Refraction  Angle  (deg) 


Range  (m) 


I  t2  /i  t2 

Figure  15.  Normal  mode  energy  fluctuation,  |4,|  /lAol  >  angle  of  45°  and 

mode  2,  overlaid  with  the  refraction  angle  in  degrees.  The  75  Hz  case  is  shown  in  the  left 
panel  and  150  Hz  on  the  right. 
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Figure  16.  Normal  mode  energy  fluctuation,  |4,|  /|4,o|  >  75  Hz  for  mode  2  and 

ISW  angles  of  5°,  15°,  35°  and  65°.  The  fluctuation  pattern  for  the  5°  ISW  angle  is  shown 
in  the  upper  left  panel,  15°  upper  right,  35°  lower  left,  and  65°  lower  right. 
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150  Hz  Mode  2  Energy  Fluctuation  5°  ISW  Angle 
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150  Hz  Mode  2  Energy  Fluctuation  35°  ISW  Angle  (|a^p/|a^p) 
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Figure  17.  Normal  mode  energy  fluctuation,  |4,|  /|4,o|  >  150  Hz  for  mode  2  and 

ISW  angles  of  5”,  15°,  35°  and  65°.  The  fluctuation  pattern  for  the  5°  ISW  angle  is  shown 
in  the  upper  left  panel,  15°  upper  right,  35°  lower  left,  and  65°  lower  right. 
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75  Hz  Mode  2  Energy  Fluctuation  !SW  Grazing  Angle  vs  Cross  Range  x=2500m  75  Hz  Mode  2  Energy  Fluctuation  ISW  Grazing  Angle  vs  Cross  Range  x=3500m 


ISW  Grazing  Angle  (deg) 

Figure  18.  Normal  mode  energy  fluctuation  as  a  function  of  ISW  angle,  , 

at  75  Hz  for  mode  2  and  ranges  of  2.5,  3.5,  and  4.5  km.  The  fluctuation  pattern  for  the  2.5 
km  range  is  shown  in  the  upper  left  panel,  3.5  km  upper  right,  and  4.5  km  lower  left. 
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150  Hz  Mode  2  Energy  Fluctuation  ISW  Grazing  Angle  vs  Cross  Range  x=2500m  150  Hz  Mode  2  Energy  Fluctuation  ISW  Grazing  Angle  vs  Cross  Range  x=3500m 


ISW  Grazing  Angle  (deg) 


Figure  19.  Normal  mode  energy  fluetuation  as  a  funetion  of  ISW  angle,  , 

at  150  Hz  for  mode  2  and  ranges  of  2.5,  3.5,  and  4.5  km.  The  fluctuation  pattern  for  the 
2.5  km  range  is  shown  in  the  upper  left  panel,  3.5  km  upper  right,  and  4.5  km  lower 
left75  Hz  /  65°  ISW  angle  /  mode  2. 
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